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Abstract. In scattering theory the far field pattern describes the directional 
dependence of a time-harmonic wave scattered by an obstacle or inhomoge- 
neous medium, when observed sufficiently far away from these objects. Consid- 
ering plane wave excitations, the far field pattern can be written as a function 
of two variables, namely the direction of propagation of the incident plane 
wave and the observation direction, and it is well-known to be separately real 
analytic with respect to each of them. We show that the far field pattern is in 
fact a jointly real analytic function of these two variables. 
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1. Introduction 

One of the main themes of inverse scattering theory for time-harmonic acoustic 
and electromagnetic waves is to recover information on obstacles or inhomogeneous 
media from knowledge of far field patterns of scattered fields excited by plane wave 
incident fields. Sufficiently far away from the scatterers every scattered field has 
the asymptotic behavior of an outgoing spherical wave, and the far field pattern 
describes its directional dependence. Taking into account the incident field, the 
far field pattern u°°{x;6), x,6 G S'^""'^, where S'^~^ denotes the unit sphere in 
W^, d = 2,3, is a function of two variables — the observation direction x and the 
direction of propagation of the incident plane wave. This function is well known 
to be separately real analytic, i.e., for each 9 g S'^"^ the function Uoo{ ■ \ 9) is real 
analytic on 5^^^^, while for each x € S'^~^ the function Uoo{x; ■ ) is real analytic on 
S'''~^ (see, e.g., Colton and Kress [F). In this work we establish that the far field 
pattern is in fact jointly real analytic on S"^"^ x S'^~^. 

To this end, we consider the special case of scattering from an inhomogeneous 
medium (see, e.g., Colton and Kress [4"), but our approach immediately extends to 
obstacle scattering problems as well. The main idea is to utilize a factorization of 
the so-called far field operator, which maps superpositions of plane wave incident 
fields to the corresponding far field patterns, to show that the far field pattern 
can locally be extended to a separately holomorphic function. Then the assertion 
follows from Hartogs' theorem [S], which is a fundamental result in the theory of 
several complex variables. A similar approach has previously been employed in a 
related two-dimensional setting of electrical impedance tomography |7j . 

The joint real analyticity of the far field pattern has an interesting consequence 
for the classical uniqueness results for inverse medium scattering problems due 
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to Bukhgeim [2], Nachman [Jjj, Novikov [T4|, and Ramm [15], which state that 
(under appropriate regularity assumptions) the index of refraction of a compactly 
supported inhomogeneous object is uniquely determined by the knowledge of the far 
field pattern Uoo{x;6) for all {x,9) G S'^~'^ x S'^~^. The separate real analyticity 
implies that it is actually sufficient to know the far field pattern on some open 
subset of 5'''"^ X S'''"^ (see, e.g., f^). From the joint real analyticity shown in this 
work it follows that the far field pattern on S"^"^ x S'^~^, and thus the index of 
refraction, is completely determined by only knowing Uoo and all of its derivatives 
at a single point {x,6) e S"*"^ x 5"*"^ 

The outline of this article is as follows. In the next section we introduce our 
mathematical setting and state the main result, which is then proved in Section [S] 

2. Mathematical setting and main result 

We consider scattering of time harmonic acoustic waves by an inhomogeneous 
medium in W^, d — 2,3, caused by a plane wave incident field u^{x;9) :— e''^^'^, 
X £ M'^, with the incident direction 9 £ S"*~^, satisfying the homogeneous Helmholtz 
equation 

Au\-;9) + k^u\-;9) = in M'* , 

where fc > is a (fixed) wave number. The scattering properties of the medium are 
characterized by its refractive index n £ L°°{R'^), and we assume that Re (n) > 
and Im (n) > almost everywhere in R'* as well as that supp(n — 1) C K'' is 
compact. Then the direct scattering problem is to determine the total field u{ - ;9) ~ 
u'{-;9)+u''{-;9) £ -^^^(M'^) satisfying 

(1) Au( ■■,9) + k^nu{ ■■,9) = in M'* , 

such that the scattered field u'^{ - ;9) satisfies the Sommerfeld radiation condition 

(2) hm rf'^-i)/^ ( - iku%- ;9)) = , 

r— >oo \ or ' 

where r := \x\ and holds uniformly with respect to x := xjr £ S^~^ . It is 
well known that ([I])-© is uniquely solvable (cf., e.g., [13 Sec. 6.2]). In fact, the 
unique weak solution belongs to iJjQ^(R'') and is smooth away from supp(n — 1) 
due to standard interior regularity results for elliptic partial differential equations 
(see, e.g., 12, p. 125, Thm. 3.2]). 

The scattered field has the asymptotic behavior 

uniformly with respect to i G S'^^^ , where the jar field •pattern Uoo{ ■ \ 9) is given 

by 

(3) u^{x- 9) = cdk^ I (n(y) - l)e-''=*-^u(y; 9) Ay, x£ S"-^ , 

with Cd = e'^/V\/8^ if d = 2 and Cd = l/(47r) if d = 3 (cf., e.g., [H Sec. 6.2]). 
Since supp(n — 1) is compact by assumption, the domain of integration in ([3]) is 
actually bounded. 

We are interested in analyticity properties of Uoo as a function on S''*"^ x 5^*"^. 
To this end, we first recall the definition of a real analytic function on a real analytic 
manifold (see, e.g., Krantz [11]): A subset C M'' is called a manifold of dimension 
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m < d if for each x <E M there exists a neighborhood U — Ux oi x in Ai, and an 
open subset W C M™, and a homeomorphism (f) : U ^ 4>{U) — W. The pairs {U, 0) 
are cahed charts, and an atlas is a family {{Ua, 4>a)}aeA for some index set A such 
that {C/ajae^ is an open covering of A4. The manifold is said to be real analytic, 
if all transition maps 

are real analytic. Accordingly, a mapping / : — > C is called real analytic if 
f o (f>a'^ is real analytic for all a G A. 

Of course, a real analytic function / : — )• C is completely determined by 
its Taylor coefficients with respect to the local coordinates and therefore by the 
derivatives 

{{D^for')){y)\veK} 

for an arbitrary {U,(f)) £ {{Ua,4>a)}aGA and y G ipiU). 

Example 2.1. The unit sphere S'^~^, d — 2,3, and therefore also S"^^^ x 3"^^^, is a 
real analytic manifold: An atlas consisting of two charts is given by 

U+ = {xe S"*-! I Xd < 3/5} , U^={xe S'^-^ I Xd > -3/5} , 

and accordingly 

(a;i, ...,xd)^ ( 1 , 

\l - Xd 
\ L + : 



^ ( i r'iu,ii2 '--- MTiu,ii2 ' !l,!!2 , ^ ) ' y = {yu...,yd-i) e 



xrf 1 + 

It follows immediately that (j)+{U+) = (t)^{U^) = i?2(0), where i?2(0) C 
denotes the open ball of radius 2 around the origin, and that 

2yi 2yd-i \\y\?~l 

+ l + ||j/||2' ||y||2 + l 

,-1/ N / 2yi 2yd-i l-\\yf\ , \ ^ m>d-i 

\l + \\y\\^ l + \\yV BII +I^ 
Since the transition map 

(0_ o ct^-^){y) = ^ ^ (0+ o ct^-J){y) , y G 52(0) \ {0} C R""-' , 

\\y\\ 

is real analytic and S'^^^ = U+ U we have equipped S*"^^^ with a real analytic 
structure. 

For later reference we note that canonically extends to a bounded function 

: y — > that is separately holomorphic with respect to each variable zj, 
j — I, . . . ,d — 1, with V C C''^"'^ being a sufRciently small open neighborhood of 
^2(0) C R'^^i in C''"^ (Here the map K'^-i Syi-^llypgRis extended as 
9 2 H> ^ 2;| e C, not as the squared norm of C''^"'^.) o 

From the representation formula it follows immediately that for each G S'^'^ 
the function Uoo{ ■ ;0) is real analytic on S"^^^. By reciprocity, 

Uoo{x]6) Uoo{-0;-x) ioT x,9 e S"^^^ 

(see, e.g., [TUl Sec. 6.3]), and thus for each x S 5"*^^ the function ■Uoo(x; • ) is real 
analytic on 5'''"^ as well. However, the joint real analyticity of Uqo with respect 
to both variables, i.e., that Uoo locally coincides with its Taylor series in the local 



4 



ROLAND GRIESMAIER, NUUTTI HYVONEN, AND OTTO SEISKARI 



coordinates corresponding to {x, 9), is a stronger result: a standard counter example 
is the function (cf., e.g., jU p. 27]) 

x'^ + 

which is separately real analytic on M^, but not even continuous at the origin. 
To the best of our knowledge, the joint real analyticity of Moo with respect both 
variables has not been reported in the literature so far. As our main result, we 
establish this property of the far field pattern in the following theorem, which will 
be proved in Section |31 

Theorem 2.2. The far field pattern Uoo introduced in ^ is a (jointly) real analytic 
function on S'^^^ x S"^"^. 

In particular. Theorem 12.21 implies that any uniqueness result for an inverse 
scattering problem related to (Il])-([2]) requiring the knowledge of Uqo everywhere 
on 5"*"^ X S'^~^ remains valid if the values of Mqo and of all of its derivatives are 
available at a single point {x,9) E S"^^^ x S"^^^. 

Remark 2.3. It follows from a result by F. E. Browder TJ Thm. 2] that the infinite 
differentiability of Moo on S'^~'^ x S'^~^, which is more straightforward to establish 
than the joint real analyticity, guarantees that Moo is real analytic on some open 
dense subset of S***^^ x S''*~^. However, this does not ensure that Mqo is jointly 
real analytic, say, along the subset {{-0,0) \ e S"*"^ of S"*"^ x S"*~\ which 
corresponds to the important case of backscattering. 

One may also argue that the application of [T} Thm. 1] would make the proof 
of the joint analyticity in the following section even shorter. However, we feel that 
a proof based on Hartogs' theorem and a standard factorization of the far field 
operator is a more fundamental approach. O 

3. Proof of the main result 

One of the main ingredients for the proof of Theorem l2.2l is a suitable factoriza- 
tion of the so-called far field operator F : L^{S'^-'^) L^{S'^-^), 

(4) {Fg){x) / Uoo{x;0)g{0) AS{e) , 

where dS" denotes the usual surface measure on S'^~^. Since the kernel of this in- 
tegral operator is bounded it follows immediately that F is compact. Furthermore, 
F can be decomposed as a product of three simpler operators as outlined in the 
following lemma (see, e.g., [S] or [31 Sec. 4.3] for related factorizations). 

Before stating this result, we introduce the three operators appearing in the fac- 
torization. To this end let i? C M'^ be a bounded domain such that supp(?i — 1) C B. 
Then the Herglotz operator A : L'^{S'^^^) — > L?[B) is defined by 

(5) {A^){y) / cl5(0) , yEB, 

and the corresponding adjoint operator A* : L^(B) — > L'^{S'^~^) is given by 

(6) iA*c^){x)= [ ^iy)e-"'^-ydy, x € S^-' . 
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Moreover, given / G L^{B) we consider the source problem 

(7) Aw + fc^™ = -fc2(n- 1)/ inR''. 

Here and in what follows, we interpret L^{B) as a subspace of L'^{W^) by identifying 
the elements of L'^{B) with their zero continuations. As in [lOl Sec. 6.2], it can be 
seen that ([7]) together with the Sommerfeld radiation condition ([2]) has a unique 
weak solution v G H^^^CM.'^). Moreover, it easily follows that the linear operator 
T : L^{B) -> L^{B), 

(8) Tf -.^k^n -!){} + v\b), 

is bounded and the far field pattern v^o & L'^{S'^^^) of v depends continuously on 

feL^iB). 

Lemma 3.1. The far field operator F can be factored as 

F = CdA*TA 

with A, A* , and T as in and ([8]), respectively. 

Proof. To begin with, we define a bounded operator W : L^{B) — )■ L'^{S'^^'^), 

Wf := Woo , 

where Wqo is the far field pattern of the solution to ([7]) with © . Observe that any 
solution u — + u'^ oi ^ satisfies 

Au'i ■■,0) + k^nu'i • ; 6*) = -fc^(n - l)u\ ■ ; 0) in W^, 

and thus it follows via superposition from the definition oi F in ^ and of A in ([5]) 
that F = WA. 

For any 4> e L^iB), it holds that A*4' = c^^Woo, where w G Hl^^{R'^) is the 
unique weak solution of 

Aw + k'^w = -(j) in R'^ 

coupled with In particular, if = Tf e L^iB) for some / G L'^{B), then 

Aw + k^w = -Tf = -k^{n - 1)(/ + v\b) = Av + k^v in R'^ , 

where v is still the solution of ([7]) with l^. Consequently, u = w in R'', and we 
conclude that A*Tf = c^^Voo = c^^Wf for any / G L'^{B). Altogether we have 
established that F = WA = CdA*TA, which completes the proof. □ 

Since A is an integral operator with a smooth kernel, it can be extended to 
a bounded operator from H~''{S'^~^) to L'^{B) for any s > 0. Similarly, A* is 
bounded from L'^{B) to iJ'*(5'^-^), s > 0. Due to the denseness of L'^{S'^~'^) in 
H~''{S'^~^) for any s > (cf., e.g., [12])j A* remains the adjoint of the extended A 
in the sense that 

(A^)^dy = {^P,A^)^,_„ ^ G H-'iS"-^), ^ G L2(B), 

where (•, •) gd-i is the (bilinear) dual evaluation between {S'^~^) and H''{S'^~^). 
In particular, F extends to an operator from Si'{S'^~^) to 2i{S'^~^) that is bounded 
from R-'iS'^-'^) to i/"(S"*-i) for any s G M. 
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Corollary 3.2. The Jar field pattern Uoo can he represented as 

J B 

for all {x;e) € S'^-^ x S'^-\ 

Proof. Denoting by 6^ S 7i-('i-i)/2-e(5d-i)^ e > 0, the delta distribution with 
singularity at z on S'^~^, we may write 



Uooix;9) = {Sx,FSe)sd-i ^ Cd \ {A5x){TA50)Ay = Cd 
and the proof is complete. 



B 



-ikx-y I 



□ 



Now, let ([/, 0) be a chart in the analytic atlas of 5** ^ as described in Exam- 
ple O We define g : (f>{U) x K'' ^ C by 

g{x,y)^e-"^^-'^^yy. 



Using the holomorphic extension 



y — >■ of ^ from Example 12.11 we 



immediately obtain the corresponding extension g : V x 



of g. Obviously, g{z, y) is holomorphic with respect to z G V for any y £ M."^. 

Lemma 3.3. Let G : F — > L^{B) he given hy G{z) — g{z,-). Then G is holo- 
morphic separately in each variahle Zj, j = 1, . . . ,d — 1, when the other variahles 
assume arhitrary fixed values. 

Proof. Pick an arbitrary 1 < j < d ~ I and let M > be such that 



< M 



-<d-l 



and y E B. 



for all z = (zi, . . . , Zd-i) in the bounded set V C 

Fix z G V and let r > be such that {w G C^^^ : \w — z\ < r} C V . Denoting 
for any 77 = (0, . . . , 0, ?7j, 0, . . . , 0) G C''"^ by F C C*"^ the straight line between z 
and z + rj, we find that 



g{z + rt,y)-g{z,y) dg 



Vj 







< 


-1 




V] Ji 



dl 

dzi 



{z,y))dx <-M\vj\ 



for all y e i? and 7^ 77^ G C such that |77j| < r. Thus, 



G(z + 77)-G(z) dG 



dz. 



(^) 



< -M|?7,VI-B| , 



which means that G is holomorphic in Zj. 

Finally, we provide a proof for our main result. 



□ 



Proof of Theorem ] 2. 2i Suppose x,9 E S"^ ^ and let {Ux,4>x), {Ug,4>g) be charts 
as in Example 12.11 such that x G Ux and 6 E Ug. Considering the holomorphic 



extensions 



-.V ^ C*, we define Gx,G'g : V L^{B) by 
e-"='^^'(^)-«, [G'g{w)]{y) = {[Ge{w)]iy))-' = .^'^^e'^-v 



[Gx{z)]{y)^e [^e[wmy) = [[^e[wj\[y)) ^ = e- 

Due to Lemma [3.31 Gx and G'g are holomorphic separately in each variable when 
the other variables assume arbitrary fixed values. 
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As T : L^{B) — > L^{B) is continuous, the bilinear form 

{p,q) := / pTqdy 
J B 

is bounded from L?{B) x L^{B) to C, and it easily follows that the (local) extension 

JB 

of the far field pattern is holomorphic vaV xV with respect to each of the coordi- 
nates separately. By Hartogs' theorem (see, e.g., [TTl Thm. 1.2.5] or [SI Thm. 2.2.8]), 
it is therefore a (jointly) holomorphic function vaV xV , and thus it is also (jointly) 
real analytic on [V xV)C^ {W^-'^ x W^-^) (cf. [HI Cor. 2.3.7]). 
Consequently, the local representation 

of the far field pattern is a (jointly) real analytic function on -62(0) x -82(0) C 
M.'^~^ x R"*"^. Since x,9 G S'^~^ were arbitrary, this shows that Uoo is (jointly) real 
analytic on S'''^ x 5"^-i. □ 
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